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Abstract 

We study the long time motion of fast particles moving through 
time-dependent random force fields with correlations that decay rapidly 
in space, but not necessarily in time. The time dependence of the aver- 
aged kinetic energy (^p^ (i)) /2 and mean-squared displacement (^q^ (t)) 
is shown to exhibit a large degree of universality; it depends only on 
whether the force is, or is not, a gradient vector field. When it is, 
(p^ (t)) ^ i^/^ independently of the details of the potential and of the 
space dimension. Motion is then superballistic in one dimension, with 
(g^ (t)) ^ t^'^^^, and ballistic in higher dimensions, with (^q^ (t)) ~ i^. 
These predictions are supported by numerical results in one and two 
dimensions. For force fields not obtained from a potential field, the 
power laws are different: (p^ (t)) ~ t'^^^ and ((7^ (t)) ^ t^/'^ in all 
dimensions d > 1. 
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1 Introduction 



We study in this paper the motion 



q{t)=F{q{t),t) 



(1.1) 



of fast particles in random force fields with correlations that are short-range 
in space, but not necessarily in time. We consider models of two different 
general classes. In the first, upon which we focus most of our attention, the 
force is assumed to be of the form 



where the /at are smooth functions of compact support in a ball of radius 1/2 
centered at 0, with additional characteristics detailed in Section [JJ i,a > 
are a length and a time scale. The /at model a random or periodic array of 
identical, randomly-oriented scatterers, centered at points qjy, that evolve 
periodically or quasi-periodically in time. We assume ini^^M H^'at — ^Af|| > ^ 
so that the local forces /at ( , ^) do not overlap. As a result, the particle 
interacts with at most one scatterer at a time, and otherwise travels freely 
between collisions. The model therefore describes an inelastic "soft" Lorentz 
gas, i.e., a distribution of soft scatterers centered at the points qn, off which 
the particle bounces inelastically. We introduce randomness in the initial 
data, and assume the system to have finite horizon, so any trajectory of a 
free particle intersects the support of F at some future time t. 

In the standard Lorentz gas, scattering is elastic, and scatterers are iden- 
tical hard unchanging obstacles centered at fixed points q^, with a spatial 
distribution chosen either randomly, periodically, or quasi-periodically. Un- 
like the current model, energy in the Lorentz gas is conserved and particle 
motion is diffusive |BSC90] . Another diffusive model, related both to the 
Lorentz gas and to those considered here, was studied in [SPB061 ILBPj : 
there the scattering mechanism was provided by a one-dimensional peri- 
odic array of oscillators representing environmental degrees of freedom of 
the medium. The Hamiltonian interaction of the particle with the oscillator 
bath then provides, in addition to a random force, an effective friction force 
that allows it to dissipate the energy it gains, and to thus equilibrate with 
its environment. 

The force (II. 2p considered in the present paper can be obtained from 
those of [SPBOGj ILBPj by switching off the friction component of the force 
provided by the particle's back reaction with the medium. The stochastic 
acceleration exerted on the particle by the random force field leads, then, to 
an unbounded acceleration of the particle, and it is of interest to compute 
the power laws associated, e.g., with the growth in time of the particle's 
average kinetic energy <^p^ (t)) /2 and mean-squared displacement (g^ (t)). 




(1.2) 
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In the other class of models that we consider, the force F (q, t) is modeled 
as a space and time homogeneous random field satisfying 

{F {q, t))=0, {F {q, t) F {q' ,t')) = ^^C , ^) , (1.3) 

where C is a matrix function of rapid decay in the spatial variable, but 
need not decay in the time variable. For these models, as with (|1.2|) . we are 
interested in characterizing the asymptotic growth of (i)) and (g^ (t)). 

There has been a fair amount of work reported in the physics and math- 
ematical physics literature on problems of this type, partially motivated 
by questions in plasma physics, astronomy, and solid state physics (see for 
example |Stu65t ISG731 |Eij97[ IBMV08| ) . Previous mathematically rigorous 
work has mostly dealt with deriving, under suitable scalings, Fokker-Planck 
equations for the particle density (as in [PV031 iGROSj ). Unfortunately, 
analyses of this type do not directly give information about the asymptotic 
behavior of the particles' kinetic energy or mean squared displacement. The 
theoretical physics literature is mostly concerned with Gaussian random 
potentials and contradictory claims have been made regarding the power 
law growth of (p^ (t)) and (^q'^ (t))- For potential fields that are delta cor- 
related in time, but not in space, it is generally agreed (see for example 
[JK82] ) that in the weak coupling hmit (p^ (t)) ~ t and (g^ (t)) ~ t^, but 
there is some controversy on what happens when the Gaussian potential 
field has temporal correlations of nonzero and finite duration. For this case 
it is argued in [GFZOTl ILMF951 IKos92j that in d = 1, (p^ [t)) ~ t^/s 
that (^q'^ (t)) ~ t^^/^ (compatible with numerical and theoretical results pre- 
sented here). In |Hei92j . on the other hand, it is claimed that for d = 1, 
{q^ (t)) ~ t^, as in the case when the random potentials are delta corre- 
lated in time. For d > 1 it is found in |GFZ91] that {p^ (t)) ~ t^/"^, and 
that (g2(t)) ^9/4. In |Kos92| . the conclusions of |GFZ91j for d > 1 are 
contested and it is argued that for Gaussian random potentials with fast 
decaying spatial and temporal correlations (p^ (t)) ~ t^/^ in all dimensions, 
and (g2 (t)) r^t^ ioTd> 1. 

Although there is some numerical work |LMF95] that supports the pre- 
dictions for d = 1 of [( ;FZ91[ ILMF951 [Ros92] . to the best of our knowledge 
no numerical simulations have been performed in higher dimensions. To 
help resolve the existing controversy on this subject we present in this pa- 
per numerical results in one and two dimensions on a particularly simple 
(non-Gaussian) model whose random force can be expressed as in (|1.2p . and 
which allows for an efficient numerical integration of the equations of motion 
out to very long times. Full details of the numerical calculations are pre- 
sented in Section [71 but our essential results for the case in which the force 
F is derived from a potential field are presented in Figures [T] and [21 where 
we plot the quantities (v^) = ({pa/tf'i and (y^) = ({q/£)'^\, as functions 



3 




Figure 1: Numerically determined values of {y'^ (r)) and (u^) in one di- 
mension (top) and for a two-dimensional hexagonal lattice (bottom), for the 
model described in Section [71 The different symbols correspond to different 
initial conditions, as indicated. 

both of the dimensionless time t = t/a and of the collision number n, which 
labels the number of scattering centers visited by the particle. 

Our numerical results indicate that in both one and two dimensions 

{v' {r)) r^l^ {vl)^n^l\ (1.4) 

which is in agreement with [GFZ911 ILMF951 IRos92j . In one dimension the 
particle's mean-squared displacement is superballistic, with 

(/(r)>~rl2/^ {yl)^n\ (1.5) 

In two dimensions, however, {^y^ (r)) becomes ballistic, i.e., 

{y'{r))^r\ {yl)^n'/\ (1.6) 

This is different from what was predicted in [GFZDlj for Gaussian potentials, 
but in agreement with predictions made for this case in |Ros92j . 

To understand our numerical results in one and two dimensions, and to 
more firmly establish what happens for the models of the type (jl.2p and 
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Figure 2: Numerically determined values of (y^ (r)) and (y^), in one di- 
mension (top), and for a two-dimensional hexagonal lattice (bottom), for 
the model described in Section [71 In each plot, different symbols correspond 
to different initial conditions, as indicated. 

()1.3p in higher dimensions, we present in the bulk of this paper a unified 
mathematical analysis that captures the essential physics of the problem. It 
provides in particular a means for calculating the power law growth of the 
mean kinetic energy and the mean-squared displacement associated with an 
ensemble of particles moving in time-dependent random force fields of the 
types described above. 

The analysis is based on consideration of the typical trajectory of a par- 
ticle moving in a fluctuating force fleld described by (|1.2p or (jl.3p . which 
can be viewed as a sequence of isolated scattering events. We argue, in fact, 
that the motion is well approximated by a coupled discrete-time random 
walk for the particle's momentum and position. Each time step corresponds 
to one collision of the particle with a single scatterer, or to one traversal by 
the particle of a distance of the order of the correlation length of the po- 
tential. Momentum increments are treated as independent random events 
whose magnitude depends upon the particle speed. Theoretical analysis of 
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the resulting random walk reveals that the high velocity behavior of the mo- 
mentum change of the particle during one such scattering event completely 
determines the asymptotic properties of the motion. As we show, this high 
energy behavior is insensitive to the details of the force field, notably to its 
statistical properties or to the precise geometry of the scattering centers; 
the asymptotic behavior of the motion is therefore quite universal, and in 
particular not a result that arises only with Gaussian potential fields. 

Indeed, for general force fields obtainable as the gradient of a potential 
field, we find (Theorem 14. 1|) that the energy change incurred by a particle of 
velocity v satisfies AE ~ ||^^|| This fact, combined with our analysis of the 
resulting random walk in momentum and position space leads to an increase 
of (t)) that is in all dimensions of the form observed in Figure [U and as 
described by (jl.4p . In one dimension, (^q^ (i)) is predicted by our analysis to 
grow in time as in (II. 5|) , and as observed in the top left panel of Figure [2j 
In all higher dimensions it is predicted to grow as observed in the bottom 
left panel in Figure O and as described by ()1.6p . This slower growth of 

(t)) higher dimensions arises from the fact that the particle can now 
turn while traveling, as its velocity vector performs an orientational random 
walk resulting from small random deflections. 

Our analysis can also be applied to the case where F{q, t) does not derive 
from a potential field, a situation which has attracted some attention in the 
mathematics literature. We find for a non-gradient force field that the en- 
ergy change in a single scattering is considerably larger than in the gradient 
case : AE ~ 1 (Theorem 15. ip . Consequently, we predict a larger rate of ac- 
celeration (p^ (t) ) ~ t^/^ (see (15. 6p ) , that confirms rigorous results that have 
been obtained for d > 4 in |DK09l IKP79j under suitable technical conditions 
on the forces fjy in ()1.2p . Our analysis then leads to the prediction that in 
all dimensions particle motion in the presence of a non-gradient random 
force field is superballistic with (i)^ ~ t^^^ (see (|5.7p ). The fundamen- 
tal reason for the difference with the gradient field case is that the particle 
turns more slowly while traveling, because it accelerates more quickly, and 
so is less easily deflected. The difference between the two situations can be 
traced to the fact that time-dependent gradient force fields produce smaller 
changes in the particle's energy than non-gradient force fields do. This is a 
remnant of the energy conservation that is a characteristic feature of time- 
independent gradient fields. 

The rest of the paper is organized as follows. In Section [2] we introduce 
a random walk description of the motion of a particle moving in a field of 
scatterers. General features of the walk that pertain to both gradient and 
non-gradient force fields are derived in Section O Section U] is devoted to a 
derivation of the above power laws for the case of a gradient force field, and 
Section analyzes the non-gradient case. In Sectional we adapt our analysis 
of Sections [3]l5] to random force fields as described by (II. 3p . obtaining results 
for the gradient and nongradient case identical to those found, respectively, 
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in Sections [H and O Details of our numerical calculations, the results of 
which are presented in figures distributed throughout the paper, are given in 
Section [71 Proofs of mathematical results used for the analysis in Sections [3]- 
[6] comprise the Appendix. 

2 Particle in a field of scatters: a random walk 
description 

We first describe precise conditions on the functions /at in (|1.2|) under which 
we work. We systematically use rescaled variables {i > 0,a > 0) 

r = ^eM, y(r) = ^eM^ y (r) = y (r) = jp (t) , xn = ^ 

and suppose fjy to be of the form 

In (y, r) = -^cnMnQ (M'^y, ujt + (t>%) . (2.1) 

The locations xn,N G Z'^ of the scattering centers can be chosen either 
randomly (with uniform density) or lying on a regular lattice. The coupling 
constants cm are independent random variables taking values in [—1, 1] and 
distributed according to a common probability measure v not concentrated 
on 0. The Mjy are rotations belonging to SO((i, M) and are also i.i.d., ac- 
cording to the left-invariant Haar measure on SO (d, R). Thus, the scatter- 
ers are identical objects randomly oriented in space, all described by the 
same function g : x T™' — > R'^ which is smooth and supported in the 
ball of radius 1/2 in its first variable; T™ = R'^/Z™ is the m-torus and 
u G R™, 1 1 a; 1 1 = 1. When u) has components that are independent over the 
rationals, the force is quasi-periodic in time, otherwise it is periodic. The 
parameters (fPj^ G are i.i.d. random initial phases, uniformly distributed 
on the torus. We write d/U (M, c) for the above described probability mea- 
sure on SO (d, R) X T"^ x [—1,1]. The force may or may not derive from a 
potential. The above class of models is sufficiently rich to allow for the de- 
scription of pulsing, vibrating, and rotating scattering centers; for an explicit 
example, see Section [71 

In the rescaled variables, the equations of motion (ll.l|) - (ll.2|) become 

y{T)=G {y (r) , r) , G (y, r) = ^ cnMn9 (M^' [y - xn) , ujt + 4>n) ■ 

N 

(2.2) 

One should think of g [y^i^r + (p) as the force produced by a soft, time- 
dependent scatterer centered at the origin; G then describes a field of iden- 
tical scatterers, randomly oriented, and centered at the points xn- We 
assume the system has a finite horizon, so that the distance over which a 
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Figure 3: A particle at time t„ impinging with velocity Vn and impact 
parameter 6„ on the nth scatterer, centered at the point 

particle can freely travel is less than some fixed distance L > 0, uniformly in 
time and space and independently of the direction in which it moves. Thus, 
with probability one, for all (y, such that G (y, r) = 0, 

supir > I VO < r" < r , G (y + vt" , t + t") = 0} < 

ll^ll 

We consider a particle which at time tq = is close to a scatterer at = 
and moving toward it with initial velocity vq along an initial direction that, if 
followed without deflection, would find the particle at its closest approach to 
the force center located at a point defined by the impact parameter bo G M"^ 
(see FigureE]). After inelastically scattering from the center at xq the particle 
moves freely with a new velocity vi until it encounters a second scatterer, and 
in this way it undergoes a random succession of scattering events. The nth 
scattering event begins, by definition, at time r„ when the particle arrives 
with incoming velocity Vn at the point (see Figure [3]) 

1 , , „, „ 1 

Un = Vn - -^en + bn, • 6^ = 0, ||0n||<- 

near the scattering center at yn = xp^„, where Cn = fn/H'^n 11, and the impact 
parameter bn is a vector perpendicular to the incoming velocity vector. The 
nth scatterer itself is characterized by its orientation M„ := Mn^, its phase 
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(j)n '■= iOTn + (j)j^^ at the time that the particle encounters it, and the couphng 
strength Cn ■= cn„- 

The change in velocity experienced by a sufficiently fast particle at the 
nth scattering center can be written (Proposition IA.2p 

Vn+l =Vn + R{ n ) Cn 

) (2.3) 

where, for ah w G R^, b E withw6 = 0, and (M,(/),c) E SO (d,M)xT'"x]R, 

r+oo 

R{v,b,M,<j),c) =c dT'Mg{M-^y{T') ,ujt' + (/)) (2.4) 
Jo 

in which y (r) is the unique solution of 

y{T) = cMg{M'^y{T),iOT + (l)) , y(0) = 6-J^, y{0)=v. 

^ ll^ll 

After leaving the influence of the nth scatterer, the particle then travels a 
distance r]n with velocity fn+i to scatterer n + 1, which it encounters after 
a time At„ = r/„/||z;n+i||. 

Based upon this description of the dynamics, and ignoring the role of 
recollisions, we now argue that the motion of an ensemble of particles moving 
in a force field described by ()1.2p is well approximated by a coupled discrete- 
time random walk in momentum and position space. Each step of the walk 
is associated with one scattering event, where the variables M„,(/)„,c„ that 
characterize the scatterer, and the variables r]n^ bn that characterize the ap- 
proach of the particle onto the scatterer, are drawn from distributions that 
characterize them in the actual system of interest. Thus, starting from a 
given initial condition (yo,fo), we iteratively determine the velocity, the lo- 
cation, and the time of the particle immediately before the nth scattering 
event through the relations: 

Vn+l = Vn + R{Vn,Kn) ^ 

Tn+l = r„ + pf^ I (2.5) 

Vn+l = yn + r]*en+l ) 

where = {bn,Mn,(j) mCn)- The parameters (^Mn,(pn-iCn ) are independently 
chosen from the distributions already described (the distribution for (pn be- 
ing the same as for cp'^). Without the loss of any essential physics, we have 
in ()2.5p replaced the random variable 7]n at each time step with the average 
distance rj^: = (rjn) < L between scattering events. The bn are indepen- 
dently chosen at each step uniformly from the d — 1 dimensional ball of 
radius 1/2 perpendicular to To summarize, this random walk describes 
a particle that moves freely over a distance rj^^ then meets, with random 
impact parameter, a randomly oriented scatterer at a random moment of its 
(quasi-)periodic evolution. After scattering, the process repeats itself. Our 
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basic assumption, therefore, is that this gives a good description of a typical 
trajectory in the real system. 

In what follows we write (•) for averages over all realizations of the ran- 
dom process k„. In Sections [3][6] we study the asymptotic behavior of this 
random walk, under conditions expressed in the following hypothesis: 

Hypothesis 1. g ^ C'^ (M°' x T"*) is compactly supported in the ball of 
radius 1/2 centered at the origin in the y variable. The function g and its 
partial derivatives up to order three are all bounded, and we write 

< 5max := \\g\\oD < +00. 

If g {y,4>) = —VyW {y,<j)), we suppose W ^ {W^ x T™) also is supported 
in the ball of radius 1/2 centered at the origin in the y-variable. Moreover, 
(to ■ V^) W ^ and, if d = 1, we require that, for some cj) G T™, 

/ dy {uj-V^)W{y,<t>)^Q. (2.6) 



3 Analysis of the random walk: general consider- 
ations 

We now turn to the analysis of the large n behavior of the first equation of 
(I23D 

Vn+l=Vn + R{Vn,l<'n) (3.1) 

which is independent of the others. We assume that the particles are fast, 
meaning ||vo|P >> c^max (Lemma lA.ip . For that purpose we need to under- 
stand the high momentum behavior of the momentum transfer R{vn, Hn), 
as well as of the energy transfer 

AE {Vn, Kn) = ^ (^{Vn + R (?J„, Kn)f - V^J . (3.2) 

First order perturbation theory allows one to write (see Proposition IA.2I for 
details) 

R {Vn, Kn) = ^ J dX Mng (^M'^ +(^X-^ , ^ + 

+o(IKI|-^). 

More generally, if g is sufficiently smooth, one can write, for € N, and 
{v,k) G R^'^ X SO (d,M) X T™ X M, ^--y = 0, 

= E^^ + o CM-"-'). -M- <^-=' 

k=l " " " " 
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Note that 



a 



(1) 



./;;dAM,(M-(.+(A-i).) 



and 



a 



(2) 



/+00 / 
dA \drMg ( M"^ 



6+ A 



in which we have introduced the suggestive notation 

dr ■■=uj- Vs- 



Hence 



where 



(3.4) 

(3.5) 

(3.6) 
(3.7) 





















= (««•« 


/3W 





a^^^+e-a(3)) > (3.8) 

a(2)+a(i).a{3)+e-a(5)). ^ 

It is easy to see that expansion (j3.3|) has rather different features when 
g' is a gradient vector field than when it is not. Indeed, when g = —VW, 
the first order term in the momentum transfer (j3.3p is perpendicular to the 
incoming momentum v, so that 

/3(°)(e,K) =e-a(i)(e,K) =0. (3.9) 

As a result AE ~ in that case. Moreover, one then has 

/ + 00 
d\drW{M-^{b + Xe,(l))). (3.10) 
-oo 

On the other hand, when g is not a gradient vector field, P^'^^ does not 
vanish and, as a consequence, AE ~ 1. This is the source of the different 
asymptotics for {v'^'j and {Vn) in those two cases, as we will see below. 
For later purposes, starting from (|3.ip - (|3.3p . a simple computation yields 



AE„ 



Cn + 



Rn 



+ 



{AEnY 



(3.11) 
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where i2„ = R and AE'„ = AE Hence, from ()3.8p . 



n I "n I ^ j \\ 112 ^ V " V " / / II II 3 

(3) / (3) ^ ^ ( (1) (1)"^ 

/ / 1 1 Vn II ^ ^11 "Vfi 1 1 



+o(ib„ir') + 



\\Vn\\^ 



4^)+o(iKn+o(i^) 



Here a^n^ = a*^^) (e„, k„) . We can write (5n = <^n + /^nenj (^n ' = 0, with 
(since ||en+i|| = 1 = ||e„||) 



fin = + <0 



For a function / depending on v and «; = {b,M,(l),c), b ■ v = 0, \\b\\ < 
1/2 we shah denote the average over parameters associated with a single 
scattering event as 

d6 



fiv)= I -^J dfi (M, <p, c) f {v, b, M, cP, c) , (3.12) 



where Cd is the volume of the ball of radius 1/2 in 



4 Analysis of the random walk: gradient fields 

In this section we consider the more interesting case where g = —VyW. 
The following theorem, the proof of which appears in the Appendix, will be 
essential to our results. 

Theorem 4.1. Suppose Hypothesis [J\ holds and that g = —VyW. 
(i) For all unit vectors e G M*^, 



a(i) (e) = = a(2) (e). (4.1) 

Moreover, for all v 

^E{v) = ^^+0{\\v\\-'>), [AE {v)f = ^ + O {\\v\\-') , (4.2) 



where 

B = - — -D'^ (4.3) 



d-3n 
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with 



D^ = ^ [ d<p [ dyody'o \\ yo - y'o II'"' d^W {yo, </>) d^W {y'^, </>) > 0. 

(4.4) 

In particular, for all unit vectors e E and for £ = 1,2, 3, 



(e) = 0, B = (e) and = {f3W (e)) > 0. (4.5) 

(ii) Let Vn he the random process defined by 113. 1\) and e„ = fn/lbn||- Let, 
for £ G N, = /3^^^ {en,Kn) ■ Then one has, for all n / n' G N, for all 
0<£<e <3, 



(4.6) 



Remark 4.2. (i) Note that part (i) of the Theorem does not involve the ran- 
dom walk It is a statement about the functions a^^^ (e, k) (e, k), 
viewed as random variables in k. 

(ii) The strict positivity of is equivalent to the requirement that P^^^ does 
not vanish identically. This follows from Hypothesis d and notably from 
the nonvanishing of the time derivative of the potential. This is as expected, 
since in a time-independent potential, energy is conserved to all orders, so 
certainly /3^^^ =0. In one dimension, the extra assumption i2. 6\) is needed 
to ensure / 0." indeed, when d=l, = as soon as the potential has 
a vanishing spatial average. In that case, some lower order term f3^^^ will not 
vanish and, as will be clear from the discussion which follows, this would al- 
ter the power laws of the stochastic acceleration. Such situations, which are 
easily treated using the methods of this paper, will not be considered further. 

We first establish the asymptotic behavior of (||fn|P), where Vn is the 
stochastic process defined by (|3.1|) . We start from the expansion (|3.7|) which 
yields, respectively 

lU, l|2 4 O/^W 

r'n+l 1 , ^M" I r^/w \\-7\ 



bnP ^ ll^nl 

1=1 



1=1 



\Vn+l\\ - \\Vn\ 



1=1 
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and consequently 

|2 



_^ \\Vn+l\\ _^ \\Vn+l\ 



3 HPi'^ + UPn^ 



VnW \\Vn\[ 
2 



i=l 



-4\ 



3te) + i(/3«)^ 



3/3i'^ + ^ ^ ^ + Oo {\\Vn\\-') + O dl^^nir') . 

I Vn 1 1 

(4.7) 



Here the notation Oq (||wn|| ^) means the term is O (||w„|| and of zero 
average. Introducing 



J„ = ^,.„ = =_and-, = -(^^ + -j=-(<i-2)>--. (4.8) 

we drop the error term in (j4.7p to obtain the one-dimensional random walk 
= en + ^ with (e„) = 0, (e^ > = 1 (4.9) 

in the variable £,n- Here the first term on the right hand side is the dominant 
term of zero average in (|4.7p . whereas the second term is a systematic drift 
term, and is its dominant term of non-zero average (when 7 7^ 0). 

From this simple random walk we can easily deduce the short time be- 
havior of the dynamics. Suppose ^0 >> ItI- Then, 



7 \ ^ 

= Co + n— + 2_^ek 



n-l 



Co 



fc=0 



where this approximation remains valid as long as |Cn — Col << Co- A short 
calculation shows this is guaranteed provided |1| 

n « A^* (Co) ~ Co ~ ll^'of. (4.10) 

This last relation gives an estimate of the number of collisions needed before 
the asymptotic long time behavior, as derived below, sets in. This depen- 
dence on the initial speed can be seen in the numerical results for the model 
described in Section[71 as shown in FigureHl We now turn to the asymptotic 
behavior of Cn n » (Co)- We will show that, for d > 2, and for k > —3, 

k 

~n6. (4.11) 
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Figure 4: The quantities (ll^'oll) (on the left vertical axis), and (||i'o||) 
(on the right vertical axis) as a function of ||fo|| for a particle moving in a 
hexagonal lattice (d = 2) as described in Section [3 



Note that this is indeed the behavior observed numerically for the full 
dynamics of the numerical models described in Section [71 as illustrated in 
Figured] for k = 2, and which we present in Figure [5] for k = —1 and —2. 

From a theoretical point of view, the result (j4.1ip is obvious for d = 2, 
since then 7 = and (14. 9p then just describes a simple random walk on the 
half line. More generally, looking at (j4.9p . because 7 > for d > 2, one 
certainly expects {^n) +00, as a result of the combined drift-diffusion 
implied by ()4.9p . In d = 1, 7 < and the second term then acts as a friction 
term. We will nevertheless show that for all 7 > —1/2 the friction is too 
small to alter the asymptotic behavior of (^n)- To this end we note that 

Ad = (2en + ACn) (AC„) . 

Again keeping only the dominant terms yields 

Ad = 2Un + 27 + 1. (4.12) 

For 7 > —1/2, we will now show that after rescaling the process ^ by 
n, it has a well-defined limit, which is a squared Bessel process of dimension 
6 = 27 + 1. This will establish (14. lip for those values of 7 and for all A; > —3. 
To see this, define, for s > 0, n G N, and < cr < s, 

yi") = -d, if ak = k-<a<{k + l)- = ak+i. (4.13) 
n n n 

^From this point onward we use the notation / (x) ~ g (x) to mean that there exist 
< c < C < +00 so that cf (x) < g (x) < Cf (x). 
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Figure 5: Numerical results showing the asymptotic behavior of (||fn||*^)! 
for the model described in Section [71 in one and two dimensions, with initial 
conditions and k values as indicated. 

Multiplying (I4.12p with s/n one finds that 

71—1 , 

yi") = + 2 Y: ^jY^^ABi-^ + (27 + 1) 

fc=0 

where 

I — k-l 

(n) 

Taking the limit n — > +oo and writing Yg = lim„__^+oo Yg , one finds 

Y, = Yq + 2 [ ^/YsdB, + (27 + 1) s, 
Jo 

where Bs is a one-dimensional Brownian motion since the e„ are i.i.d. In 
other words, the limiting process Yg satisfies the stochastic differential equa- 
tion 

dYs = 2^/YsdBs + (27 + 1) ds, (4.14) 

of the squared Bessel process of dimension 6 = 2^ + 1 (see |RY99| . Chap- 
ter 11), and is therefore a squared Bessel process. 
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Figure 6: Behavior of the cohision time {Tn) for the model of Section [3 both 
in one and two dimensions, with initial speeds as indicated. 

Thus, since ^^/n converges, we can approximate its distribution by that 
of Yi and conclude that, for all I > —1, ~ na , which is (|4.11|) . Equation 
()4.1ip in particular yields, via ()2.5p . (see Figure [6]) 

(Tn) ~ n^/' (4.15) 

and, finally 

(r;2(r)>~r2/5, r » r, (||t;o||) := ^^j^ ~ ||z;of ■ (4.16) 

We note that the asymptotic behavior does not depend on the initial 
speed ||fo||; the time scale r* (||t'o||) = ''"ATtdlvoll) which it sets in, on the 
other hand, is predicted by this analysis to grow quickly, as ||wo||^) a result 
that is verified in the 2d numerical results presented in Figure [U which shows 
the number of collisions A^* and the mean time t^, before the asymptotic 
regime is observed, as a function of ||uo||- 

For d = 1, the same power law was found for (Gaussian) random fields 
in |GFZ91j . |LMF95j . and |AWMN06j . using very different methods. For 
d > 1, the only studies we are aware of are |GFZ91j and |Ros92j . who 
deal with Gaussian random fields and who respectively find (^v"^ (r)) ~ r^^^, 
which disagrees with KWh and {v'^ (r)) ~ T which agrees with it. As 
mentioned in the Introduction, we have corroborated our predictions ()4.1ip 
and (I4.16p . including the onset of the asymptotic regime at n (Hi'oll), with 
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numerical calculations in Id and 2d, the results of which are presented in 
Figure [J and described more fully in Section [71 

Before turning back to an analysis of the full random walk (12.50 in or- 
der to determine the asymptotic behavior of the mean squared displacement 
(y^ (r)), we now first briefly discuss the validity of the assumption we im- 
plicitly made in passing from (14. 7p to ()4.9p . namely that the lower order 
terms of (|4.7|) won't alter the behavior of {v'^ (r)) that we obtained by ig- 
noring them. To get an estimate of the error made neglecting these terms, 
we will evaluate them along a typical trajectory of the random walk (|3.1|) . 
along which we showed ||f„|| ~ n^^^, and will thereby demonstrate that the 
contribution of each of the neglected terms to is smaller than n^/^, the 
contribution of the two dominant terms retained above. Note first that for 
i = 2,3 



E 



/3^*) \ 



^ A;(i-l)/6 




because (4^^) = 0, and (/^fc^/J^O = for / k' , by Theorem HHii) . 



Hence, since for i = 2,3, ^^/^ << n ~ we conclude 

(i) 

that, with the above condition on , these neglected terms do indeed 
contribute a lower order correction to ()4.1ip . The neglected term of order 
is also of zero average, and therefore treated in the same way. Unlike 
the first three terms, the error in (14. 7p that is of order ||fn||~^ = O (n^^/'^), 
need not be of zero average; however after summation over n it yields a 
contribution of order n^^^ « n}!'^ ~ and can therefore also be neglected. 
Theorem 14. 1 1 therefore implies that all neglected terms in ()4.7p provide lower 
order contributions to the asymptotics of Note the crucial role of the 

retained term in ()4.9|) involving 7, which contributes a term of exactly the 
same order as the dominant diffusive term e„ = I3n^ /D. 

We now derive the asymptotic behavior of ||y„|| and ||y (r) || (see ()4.20p 
and ()4.2ip below). We first consider the case with d > \, which clearly 
depends on how much the particle's path deviates from a straight line, i.e., 
on how much and how quickly it turns. In particular, we need to analyze the 
third equation in (12. 5p . For that purpose, we will first study the evolution of 
the unit vectors e„, which execute a random walk on the unit (d — l)-sphere. 
Note that, for all result of Theorem 14.11 ^1) and the observation that 

e„ is independent of k„, the step 5n of the walk in Cn, defined in (jS.lip . has 
a mean that satisfies {5n) = O On the other hand, the magnitude 

of the step is 

ii'^niN^+odKir^) = 11^11. 



Given that the particle has high speed at the nth collision, we now 
wish to compute how many collisions m it takes for the particle's direction 
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to change by a macroscopic amount. For that purpose, we compute the 
conditional expectation 

m— 1 m— 1 
k=0 k'=0 

We will suppose m satisfies m « (^„) ~ ||fn||^ ~ n, so that ()4.10p 
implies ||vn+m|| ~ lbn||; we will therefore approximate IliJ^+fcH by ||fri||. It 
then follows from Theorem I4.1( i) that, for all k, 

{8n+k ■ Sn+k) = ^ II ^ + O {\\Vn\r^) . 

1 1 1 1 

For the off-diagonal terms, we note that for k > k' , 

{^n+k ■ 5n+k') = {^n+k " ^n+k'^^ + l^l^n+k^in+k ' ^n+k'^ + 

In addition, the rotational invariance of the system implies that for a given 
e„+fc, the vector (^^^_jr. vanishes (see (|3.12|) for definition of the barred aver- 
age). Hence, \i k > k' 

On the other hand, writing that e„_|_fc = Cn+k'+i + ^fci rotational invariance 
also implies that the conditional expectation of given Cn+k'+i is a vector 
UkSn+k'+i of length < 2. Hence, 

(^S^+k' ■ fJ-n+k^n+k^ = (^S^+k' ■ fJ-n+ken+k'+l^ + (^S^+k' " fJ-n+kt^ken+k'+l^ , 

and 

,^n+k' ' f^n+ken+k) \ < {\^n+k' ' ^n+k'+lllfJ-n+kl 



< 3 

< - 



Sn+k' ■ [dn+k' + 5n+k'] 



jMll^n+fc'l 



o(iKir^). 



Consequently, 



\a, 



a)||2\ 



> = ^ u 114 + (Ibnir^) + m^O (||r;„||-^) . 
1 1 1 1 



Consequently, provided 

m = M, {\\vj) ~ ~ n2/3 « n (4.17) 
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Mdivji) 



Figure 7: On the left, the correlation function (e^ • eo) is plotted as a func- 
tion of m for a set of fifteen initial speeds ||i'o|| lying in the range 0.5 to 2.0. 
On the right, a numerical estimate of M^, (||uo||) obtained from the initial 
slopes of the data in the left panel are plotted as a function of ||fo||- 



we find ( 11 

Cn+m CnlP) ~ !• This shows that after iW^.dli'^ll) collisions, and 
aside from accidental cancelations between the diagonal and off-diagonal 
terms, the particle turns through a macroscopic angle with the unit vectors 
en+m covering the unit sphere. In Figure [7] we display values of M^<(||uo||) 
obtained from a numerical study of the decay of the correlation function 
(cn • Co) in the 2d numerical model described in Section [71 The observed 
power law behavior agrees with the one predicted by the random walk anal- 
ysis above. 

We now analyze the asymptotic behavior of ||yn||- For particles that 
start off with an initial speed Ht'oH, it takes typically Mi = ||vo||^ collisions 
to acquire a random direction of motion. We can then define recursively 

Mk+i=Mk + Ml, (4.18) 

from which one readily finds that ~ k^. The can be interpreted by 
remarking that, when n = M^, the particle's velocity has "turned", i.e., 
changed direction by a macroscopic amount, on average k times, while the 
trajectory along the sequence of m << M^+i — M^, collisions between 
and Mfc+i largely follows a more or less straight path. We use this picture 
to approximately compute yM^. by writing 

VMk+i = 2/Mfe + V* (^fc+i - Mk) CMf, ■ (4.19) 
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This is a rough estimate, but the idea is that, on average, the particles go 
straight for about M/^^i — steps in the direction 644 without turning. 
In view of ()4.17p and ()4.18p we can now think of these successive directions 
CMfc as randomly and independently chosen on the sphere, so that (j4.19p 
describes a random walk on a larger length scale, having independent steps 
of order rj^ (Mfc+i — Mfc) ~ A;^. This yields 

k 

Interpolating between the then allows one to write 

{WVnf) n'/\ (4.20) 

Note that, together with (j4.15p . this finally gives the result 

(||y(r)||)~r. (4.21) 

The motion of the particles is therefore ballistic in the sense that \\y (r) ||/r, 
which describes the rate at which the particle's distance from the origin 
grows, is finite on average. Note, however, that the averaged instantaneous 
speed grows as r^^^, as shown above. The particles therefore speed up, but 
turn while traveling, which decreases the rate at which they move away from 
the origin. The results of our numerical calculations for d = 2, presented 
in Figure El and described in detail in Section [7] are in agreement with the 
results of the random walk analysis outlined above. 

Finally, we briefly treat the situation for d = 1. For this case, the parti- 
cle cannot progressively change its direction, and so the analysis presented 
above does not apply. Indeed, in Id a direction change implies a complete 
reversal in its direction of motion; but this can happen only if the particle 
encounters a stretch of scatterers that causes it to completely decelerate 
first. Computations similar to the previous ones show this cannot occur on 
a time scale shorter than M^: (||vo||) ~ ll^oll^i which is the same scale on 
which, as we have shown above, the particle accelerates. On the other hand, 
for all times, we have the obvious upper bound 

(||y(r)||)< Td. {\\v{t)\\)^t'/^ (4.22) 
Jo 

It is clear, then, that for time scales over which most of the particles in the 
ensemble have not reversed direction 

(l|yn||)~n, (4.23) 

which along with (I4.15P implies (Hy (t) ||) ~ r^/^, i.e., it saturates the upper 
bound (j4.22p . At longer times, the distribution of times for the random 
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walk ()4.9p to return to the origin at ^ = 0, which governs events at which 
the velocity reverses, may alter the asymptotics. If this happens, it does so 
at times longer than we have been able to investigate numerically. Indeed, 
up to the times investigated in our numerical calculations the bound (j4.22p 
appears to accurately describe the asymptotic properties of the growth of 

5 Non-gradient force fields 

When the force field g does not derive from a potential W, we suppose the 
distribution i' of the coupling constant c is centered 

cdi/(c) = 0, (5.1) 
so that the mean force vanishes at each point y G M"^. From ()3.7p . we have 

Ibn+lf = ll^^nf + 2/3W + ^ + 1^ + O (ll^.r^) . (5.2) 

1 1 II II Vn 1 1 

Similar to Theorem 14. 1^ we now have the following Theorem, the proof of 
which also appears in the Appendix: 

Theorem 5.1. Suppose Hypothesis {1\ and i5.1\) hold. Then: 

(i) For all unit vectors e G M*^, a^'^^ (e) = = a(^) (e). Moreover, for all 

V G W^, 

A£>) = ^ + 0(||7;||-3), {vf = D'^ + 0{\\v\\~^), (5.3) 

where B' = ^'^^^^ D'^, with 

= ^ [ d</>/ dyody',\\yo - y'o\\-^'+'^ 

X (yo -y'o) ■ 9 {yo, 4>) {yo -y'o) -g (2/0, 4>) > 0. 

In particular, for all unit vectors e G and for £ = 0,1, 



/3W(e)=0, 5' = /3(2)(e) and!?' = (/3(0)(e))'>0. (5.4) 

D' > if and only if P^^^ (e, /t) does not vanish identically, which implies g 
is not a gradient vector field. 

(a) Let Vn he the random process defined by ici.l]) and Cn = Vn/\\vn\\- Let, 
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for e E N, pi^^ = /3W {en,Kn) . Then one has, for all n / n' G N, for all 

o<e<e' <1, 



Moreover, ( f/?!^^^ ) and l[in^[in^\ are independent of n. 



(5.5) 



Remark 5.2. Whether the force does or does not depend on time plays no 
role in this result, contrary to what happens in Theorem \4-l\ In other words, 
when a force field is not a gradient field, the dominant behavior of the energy 
transfer to a particle is not affected by whether it depends on time or not. 
In particular, the coefficients B' and D' do not involve time derivatives of 
the force, as do B and D. 

We now analyze the asymptotic behavior of the velocity and the position 
of the particle, as in Section^ From (j5.2p . neglecting the subdominant terms 
as in (fi71) - (fO|) . we find 

i hn\? Pn^ 1 

= + ^' ^^^^^ = = and y = - (d - 1) . 

Note that 7' > in all dimensions, so that from the analysis of ()4.9p it 
follows that (||fn||'^) ~ n'^/^. Using this in (|2.5p yields 

n 

^^«)~E7TM (IbWII)-^'/'- (5-6) 

£=0 ^ 

This is proven rigorously in |DK09j for a time- independent, non- gradient 
force field of the type ()1.2p and (I2.ip . in d > 4, under suitable additional 
technical conditions on g and the distribution of the scattering centers. 

We now show what this implies for the asymptotic behavior of ||y(T) ||. 
First, the short time scale A'^^ (^q) is now A^* (^q) ~ ~ n ~ ||fo||^- Then, 
from (13. lip we find 

(1) f (1) \ 
an - [an ■ e„ I e„ 

Cn+l =en + ^ + O {\\Vn\\-^) ■ 

\\^n\\ 

Consequently, ^||e„+m — ^nW^) ~ ~ m/n. Thus, the particle now 

turns over a macroscopic angle after M^, (||wn||) ~ ~ n collisions, 

many more than for force fields deriving from a potential (see (j4.17p and 
of the same order as the number iV* (H't'nll) ~ n of collisions it needs to 
accelerate significantly. This is simply due to the fact that the particle is 
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much faster, since ||f„|| ~ n^^'^ rather than ||f„,|| ~ n^^^, and harder to 
deflect. This reflects itself in the asymptotic behavior of ||y(T) || as follows. 
We define as before Mi = \\vo\\^, Mk+i = Mfc + M^, so that ~ 2^=, and 
VMk+i = VMk + V* (Mk+i - Mk)eMk, which integrates to (||yAfJ|) ~ M^, 
yielding 

(l|y(T)||)~r^/^ (5.7) 
independent of the dimension d of the ambient space. 



6 Homogeneous random fields 

As we now briefly indicate, the analysis of the previous sections can be 
adapted to the case where the force field is not of the form (12. 2p . but is a 
time and space homogeneous random vector field satisfying 

{G {y, r)) = 0, {G (y, r) G {y', r')) = C {y - y', r - r') . 

Note that C is a matrix-valued function, which we assume decays quickly in 
its spatial variable, but not necessarily in its temporal variable. 

In this situation, also, we expect the asymptotic motion of the particle 
to be well described by a random walk similar to the one in ()2.5p . where 
now the time step Ar^ is determined by the time the particle needs to travel 
through a distance rj^, equal to several times the correlation length (which 
equals 1 in the rescaled units used here) of the force field: 

Vn+l = Vn+ Riyn,Vn,Tn,ATn) 

T„+i = t:„ + ^, 77, > 1 I (6.1) 

Un+l = Vn+ri^en- J 

Here, R {yn, Vn, Tn, Ar„) is the momentum change experienced by a particle 
that, after arriving at y„ at time t„ with momentum Vn, travels for a time 

ATn. 

We consider first the case in which G = — VVF is a random gradient field 
such that 

{Wiy,T)) = 0, {W{y,T)W{y',T'))=K{y-y',T-T'), (6.2) 

where K isa function of compact support in B (0, 1) belonging to 
that is rotationally invariant and even in its temporal variable. 

To study the asymptotic behavior of Vn in (jO.ip we first need, as in the 
previous sections, to understand the asymptotic behavior of 

n— 1 n—1 

\Kf = Hf + Yl ^W^kf = + 2A/ffc -2iWn- Wo) , (6.3) 

A;=0 fc=0 
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where Hk = H (y^, Vk,Tk) = |bfc|P/2 + Wk, and Wk = W {yk,Tk). Introduc- 
ing 

/\H {y, V, r, At) = H {y {t + At) ,v {t + At) ,t + At) - H {y, v, r) 
we find 

AH (y, V, T, = AHi (y, v, r, +AHn (y, v, r, +0 (||?;|r^ 



where 



and 



AHi { y, V, T, ^] = ^ C dXdrW (y + r],Xe, t + ^ 



V J V 







AHji{y,v,T,^^ = j\xvdrW (^y + r],\e,T + ^^ 



We then have the same kind of result as in Theorem 14.11 the proof of which 
is immediate: 

Proposition 6.1. Under the above conditions, (a^^^) = = (a^^^) and for 
all V G M"', 

{AH (v)) = ^ + o (ii^ir') ' {i^H (v))') = - J + o (ibir 

where 

B = {d-3) - 2 (d - 4) D^ = 2 (r?,K(0) - i^^^)) , 

and 

= /' d/i {-d^K (^e, 0)) , = [' dfi {-fid^K ifie, 0)) . 
Jo Jo 

Proof. Noting that {AHj (v)) = 0, we find 

{AH (v)) = {AH I J (v)) 
= ^ £ dA J,' dX" (A - A") (AdtK) (r,, (A - A") e, ^-^^) 

= ^ /o'dA/o'dA" (A - A")' {AdfK) {r,. (A - A") e,0) + O {\\v\\- 

= ^ f,' dX (1 - A) A2 {AdfK) (r?,Ae, 0) + O [\\v\\-') . 

Using the rotational invariance of Ad^K (-,0), and integrating by parts, we 
obtain the above expression for B. □ 
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Scaling ||fn|P by (s/n)^''^ in ()6.3p and taking n to infinity, one finds that 
the limiting process satisfies the stochastic differential equation 

2dB 2( l\da 

It then follows from the Ito formula that Y^j = satisfies the stochastic 
differential equation of the square of the Bessel process |RY99] of dimension 
5 = 27 + 1 . By taking -q^, sufficiently large we can make 7 arbitrarily close to 
{d — 2) /6. The analysis of the random walk is therefore entirely analogous 
to the one in Section [U yielding in particular the same power laws for the 
growth of {v'^ (r)) and (y^ (r)) as in (fi?TH]) and (jMH)- 

In the case that G is not a gradient field we still assume it to be rota- 
tionally invariant and reflection symmetric. This implies that there exist 
functions Ai and A2 such that the correlation function is of the form 

C7(y,r)=Ai(||y||,r)Pj, + A2(||y||,r)P^, 

where Py is the orthogonal projector along the direction of the vector y and 
Py + Pj/ = Id- We in addition assume that Ai and A2 are functions 
that decay fast in their spatial variable, and that for all r G R, Ai (^r) and 
A2 (•,t) are compactly supported in [0, 1]. Under these assumption, we then 
prove an analogue of Theorem I5.lt 



Proposition 6.2. Under the conditions stated above {cS-^^^ = = (a*-^^), 
and for all v ^W^, 



{AE{v)) = —^ + O IIHr' , (iAEiv)f) = + O ^"-"-1 



where 

B' = rj,{d- 1) - (d - 2) i^'(i) , D'^ = 2 (r/,i^'(0) - i^'^^)) > 
and ^ ^ 

K'W= /" dMi(^,o), i^'W=/" d^Mi(/^,o). 

Jo Jo 

Proof. A computation of R (y, v, r, 77=,,/ \\v\\) = J^^*''/'''''! Q (r') , r') dr' to 
second order in perturbation theory gives 

R{y,v,T,r]^/ \\v\\) = Rj {y,v,T,r]^/ \\v\\) + Ru {y,v,T,r]^/ \\v\\) 

with 

i?/ (y, v,r,?7^,/ ||v||) = f-^ / dAG y + T/=^Ae,r + -f-^ 
ll^^ll Jo V h\\ 
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and 

Rii{y,v,T,7^,/\\v\\) = -Ih^ I dX I dA"(A-A")x 

\\v\\ Jo Jo 

(G {y + ??,A"e, r) • V) G (y + 7?*Ae, r) + O (||f 

Hence, (a^^'') and (a'-^^) , and consequently and (Z?'-^^) , vanish. Using 

w ■ Fy {v) = {v ■ y) {w • y) jy^ then yields 



^ /^dA [\x"{Ai (L|A-A"|,0) +(d-l)A2(r?*|A-A"|,0)) 

2 Jo JO 

7?2 /" dA (1 - A) (Ai (r/,A, 0) + (d - 1) As ir]*X, 0)) , and 







a(3) . e 

= r/2 dA (1 - A) {rj,XA[ (A, 0) + (d - 1) (Ai (A, 0) - A2 (A, 0))) 

JO 

= ??2 [ dA {{{d - 2) - A (d - 3)) Ai (7?,A, 0) - (d - 1) (1 - A) A2 (rj^X, 0)) . 

JO 

Adding the last two equations and making the change of variables ix = r]j,X 
yields the above expression for B' . □ 

Analysis of the random walk now proceeds along the lines of Section [SI 
yielding the same power laws as obtained therein. 



7 Numerical results 

To verify our theoretical analysis of the motion of a particle in random 
force fields presented in the previous sections, we performed numerical cal- 
culations for a periodic array of soft scatterers in one and two dimensions. 
For the two dimensional case we employed a hexagonal lattice, with, for 
N = {Ni,N2) et^^XN = N1U+N2V, where n = (l,0), v = \{l,V^. 
We focused on the case in which the force fields associated with the scat- 
terers were derived from a potential, taking W to be of the form of a time- 
dependent, flat circular potential, 

= /((/>)x(^M^, yeM^d = l,2, 
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where x (a^) = 1 if < x < 1 and x i^) = otherwise. Here the parameter 
satisfies ^ < y* < 1/2 to ensure the system has a finite horizon. Three 
different choices were explored for the function /, namely, 

/i (</.) = cos (27r0) , /2 (</.) = 1 + cos^ (27r0) , € [0, 1 [, 

each of which leads to a time-periodic potential, and 

/3 {<t>) = h {<l>iA2) = cos (27r^i) + cos (27r02) • 

In the latter case, the frequency vector ui was chosen to be w = (l, \/2) so 
that the resulting potential is quasi-periodic in time. The phases (f)f^ were 
chosen uniformly on the torus, independently for each scatterer. Coupling 
constants cat were either drawn independently from a uniform distribution 
on [0, 1/2], or set to a fixed value cn = 1, or cjq = —1, for all N. 

Depending on the phase and the choice of coupling constants, each such 
potential describes a centrally symmetric potential barrier or well, whose 
maximum/minimum oscillates in time. For the choice / = /i or /a, any 
given scatterer will sometimes act as a potential well, and at other times as 
a barrier, depending on the sign of c^fi {(t>N + ''") or of c^fz {4>N + ^t) at 
the time r of arrival of the particle; on average the force at a given point 
in space always vanishes. When / = /2 on the other hand, and c^v = 1 
for all CNf2 {(j) + r) is always positive, yielding a lattice of oscillating 
potential barriers, for which the average force at a given point in space does 
not vanish. Similarly, when / = /2 and cjv = — 1 for all N one obtains a 
lattice of oscillating potential wells. In all cases that we studied numerically, 
the system had finite horizon. 

Motion of a particle through an array of such scatterers can be computed 
iteratively, by using energy and angular momentum conservation at the entry 
and exit of the particle from the support of the potential, and without a 
numerical integration of a second order differential equation. This allows 
one to compute the motion of the particle numerically for very long times, 
as required to properly study the asymptotic regime. 

In our calculations, each particle was initially placed at a point randomly 
chosen on the boundary of the scatterer at the origin, with an initial velocity 
drawn with equal probability from all possible outward directions. For each 
ensemble of initial conditions, the initial speed ||uo|| of the particle was kept 
fixed (with values indicated in the figure captions, or in the figures them- 
selves). Displayed results represent averages over, typically, 10^ trajectories 
for each initial speed. For convenience of presentation, data appearing as a 
function of time r or collision number n in the figures presented throughout 
the paper represent a subset of the data generated, evaluated at values of 
the independent variable that are equally spaced on a logarithmic axis. 

A general finding of both our numerical calculations and of our theo- 
retical analysis is that the power law behaviors associated with stochastic 
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acceleration is independent of the precise form of the potential employed; 
in particular it does not depend on whether the average force vanishes or 
not. Thus, in the figures that appear in the paper we have chosen to present 
numerical results only for the case where / = /i and cat is uniformly dis- 
tributed in [0, 1/2]. 

For this specific model, Figure [U shows the evolution of the particle's 
mean kinetic energy, both as a function of time r and as a function of collision 
number n. As noted in the text, one observes excellent agreement with 
the power law behavior predicted by our analysis (see p.4|) . (j4.1ip . (j4.16p ). 
independent of dimension. One also notices in this figure that the asymptotic 
regime is reached after an initial period which ends after a number N^, (||uo||) 
of collisions that grows with ||fo||. The value of (Iboll) was computed 
numerically for fifteen values of ||fo|| between 0.5 and 2, and the result is 
presented in Figure HI The observed power law A^*(||fo||) ~ Iboll^j is as 
predicted in Section [5] (see equation (j4.10p ). 

Similarly, Figure [2] shows for the specific numerical model described 
above, the evolution of the particle's mean squared displacement as a func- 
tion of T and n. We find that the power laws obtained in one dimension 
(see ([LS]) and ([03]) ) and in two dimensions (see and ([O0]) - (fOn) ) are 
indeed different, and precisely as predicted by the analysis of Section HI 

In order to obtain analytical results for a sufficiently general class of 
potentials, the theoretical analysis of Section [J] assumed scattering potentials 
that are smooth, which the potentials used in the numerical calculations are 
clearly not. Indeed, running the numerics for sufficiently long times with a 
smooth potential would involve repeatedly solving a second order differential 
equation; this would lack precision and be too time consuming. Explicit 
computations specific to the square potential, however, show that formulas 
(|3.7p and (|3.3p remain valid, and that their dominant terms have the same 
behavior as in the analysis presented, so that our arguments go through 
unaltered. This lends additional support to our claim that it is the high 
energy behavior of the energy and momentum transfer in a single scattering 
event that determines the asymptotic behavior of the particle, and suggests 
that the results are even more universal than is implied by our analysis. 

As a closing comment we note also that in our numerical models the 
potentials are rotationally invariant, and the lattices are ordered. Thus, 
when Cat is constant, the only randomness left in the problem is in the initial 
phases (j)N the scatterers and the initial directions cq of the particles. 
Thus, the essential randomness necessary for the validity of our analysis 
arises from the dispersive nature of the scattering event itself, which leads to 
a random sequence of scattering events when evaluated along the trajectory 
that the particle follows. 
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A Proof of main theorems 



In this appendix we provide proofs of Theorems 14.11 and 15.11 We be- 
gin with some preparatory material. It is convenient to write g (y, r) = 
Mg (^M^^y,u}T + cj)), suppressing the variables (f) and M from the notation. 

The estimates below are all uniform in <j) and M. Note that when g = 
-VW, then g = -VW, with W{y,T) = W {M'^y^u^T + </)). We need to 
study the solutions of 

y [t') = eg {y (r) , r') = -cVW {y (r) , r') , y (tq) = yo, y (tq) = vq. 

(A.l) 

For any initial condition y (tq) = yo, v (tq) = wq, we define 

v± = lirn ij (r) . (A.2) 

T—*±00 

That these limits exist if ||t;o|| is large enough is a consequence of the fol- 
lowing lemma. 

Lemma A.l. Suppose Hypothesis{l\holds. LetTQ G M and suppose {yo,VQ) E 
M^*^ satisfies \\yo\\ < 1/2, ||fo|P > 12c5max- Then there exist unique Tm < 
To < Tout so that \\y (rin) || = | = \\y (Tout) ||- Moreover 

1 5 

< min{ro - Tin, Tout - To} < max{ro - Tin, Tout - To} < T. 77. (A.3) 

IpoII II^'oII 

The lemma roughly says that any particle that is at some instant tq 
inside the region where the potential does not vanish and that has enough 
kinetic energy at that moment, has entered it in the past and will leave again 
in the future, spending a time of order to cross it: both the upper and 
lower bounds in (jA.Sh will be used in the proof of Proposition IA.3I below. 
Note that the Lemma does indeed imply the existence of the limits in ()A.2p . 

Proof. From (jA.lh . 

y (t) = yo + ^^0 (t - To) + c / dr' / dr'g {y (r") , r") , (A.4) 

J TO J TO 

SO that Q {t — To) < \\y (r)|| < -P (t — tq) , where 

. (t - To)^ 1 

Q(t-To) = -Cffmax ^ h ||t^0 1| |t - To | - - 

(t-Tq)^ , 1 

IX ^ h IpoIIIt - TqI + -. 

^ - P((T-), with 




i\ J ll^oll /, 6c5, 

1 , and CT+ 



C^max \ y ll^^oll^ 
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Note that cj_ < cr_|_. Since ||y (tq ±cj_) || < | < \\y {tq ±0"+) ||, it is clear 
there exist Tin, Tout satisfying 

To - 0-+ < Tin < To - (T_ and To + fJ- < Tout < To + (T+. 

Uniqueness follows from the observation that g vanishes outside the ball 
of radius 1/2 so that the trajectory can enter and leave the ball of radius 
5/2 only once. Equation (|A.3p now follows from the observation that, if 
< a; < A-^ < 1, then 

, 1 Va , 1 Va 

V 1 + X — 1 > T", 1 — v/l — a; < 'T- 

- 2,/TTA - 2 

It is enough to choose ^4 = 1 in the first inequality and ^4 = 2 in the 
second. □ 

With v± from ()A.2p . we define, for all {yo, vo,tq) so that ||vo|P ^ 12c(7max5 
Av{vo,yo,To) =v+ -V-, AK {vo,yo,To) = ^ [v^ - v'i) . (A.5) 
Note that both Av and AK are constant along trajectories: 

'^v{vo,yo,To) = Av(t;(r ),y(r ),r'), AK{vo,yo,To) = Aif (r;(r'), ^(t'), t'). 

(A.6) 

We therefore think of them as functions on the space of all trajectories 
with sufficient kinetic energy. We are interested in understanding the high 
velocity behavior of AK and of its average over all those trajectories that 
enter the support of the potential. We will see that, when g = — VVK, 
AK ~ ||fo||~^ (Proposition lA.ST) but that the average of AK vanishes up to 
terms of order ||wo||~^ (Proposition lA.4]) . In other words, the dominant terms 
in AK vanish on average. This observation is the essence of Theorem 14. 1^ 
as we will explain below. 

Proposition A. 2. Suppose Hypothesis ll\ is satisfied. Let yo,VQ S and 
Tq £ M, with ||fo|P > 12cgniax- Then 

c /■+°° 

Av{vo,yo,To) = T. — 77 / 5(yo + Aeo,To) dA 

ll^^oll J-oo 
C f'^°° 

+]n^2 drgiyo + Xeo,To)XdX + 0{\\vo\\-^) . (A.7) 

1 1 ^0 1 1 J — oo 

The error term is uniform in yQ,TQ, c € [—1, 1] and in cq = V()/\\vq\\. 

Proof. Suppose first ||yo|| ^ ^- It follows from Lemma lA. II that, under the 
stated condition on ||uo||, there exist unique entrance and exit times Tin and 
Tout to the ball of radius 5/2, with Tout — to and tq — Tin of order ||t;o||~"^. As 
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a result, for r < Tin and r > Tout the particle executes a free motion with 
speeds v- and in the region where the force g vanishes identically. From 
(jA.ip . one now readily concludes 

/"Tout 

^v{vo,yo,To) = c g{y{T),T)dT 



HyO + Vo{T -To) ,T)dT + 0{\\vo\\ ^) , 

where we used \\y (r) — (yo + vo{t — tq)) || < §(?max (''" — ''"o)^ > which follows 
easily from (|A.4p . Let us now remark that Lemma lA.ll implies that 

\\yO + VO (^out/in - To) II > Ibollkout/in " ^ol - ^ > 1/2. 

As a result, we can extend the r integration to the full real axis; indeed, 
the integrand vanishes for r < Tin and for Tout ^ The change of variables 
^ = Iboll (t - To) yields 

^v{vo,yo,To) = T. — 77 / 5 ^0 + ^^^0,^-0 + 71 — H" d-^ + O (ll^oir^) , 
Iholl J-oo V \m\\J 

(A.8) 

so that a first order Taylor expansion yields the result. 

We now consider the case where \\yo\\ > 1/2. We may assume the particle 
trajectory intersects the ball of radius 1/2 centered at the origin: otherwise 
Av {vo,yo,To) = 0, and then the result stated certainly holds. Suppose 
therefore the trajectory intersects that ball and that yo- vq < 0. Then there 
exists a unique time > tq when the trajectory enters the above ball: so 
y (t) = yo + vo{r - To) for all t < \\y (r*) || = 1/2 and y {n) ■ vo < 0. 
Clearly 

Av {vo,yo,To) = Av {vo,y{n),n), 
and we can apply the result of the first part of the proof to write 



{vo,yo,To) 



9 i^y in) + Aeo, r, + dA + O (||?7oir^) • 

(A.9) 



The change of variables 

A = A + 11^0 II (n - To) 



transforms (|A.9p into (jA.Sp . which concludes the proof. The case where 
2/0 • Wo ^ is treated analogously. □ 

When g = —VW, we need the high ||t;o|| expansion of AK up to order 
||fo||~^ obtained in the following proposition. 
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Proposition A. 3. Suppose HypothesisUlis satisfied and suppose g = —VW. 
Then, for all vq G R"' such that ||fo|P > 12c5max o^nd for all yQ G M^, 

AK {vo, yo, To) = AKi {vo, yo, tq) + AKn {vo,yo, tq) + O (||i;oir^) (A.IO) 

where 

AKi {vo, yo, To) = ^ / dA d^W (yo + Aeo, tq + ^ ) , (A.ll) 
ll^^oll Jm V Folly 

and 

AKn {vo, yo, tq) = --nAn / dAV^^T^ (yo + Aeo, tq + 77^^ 

ll^^oll"* Jr V Foil/ 

dA'^ dX"VW (^yo + X"eo,To + ^^ . (A.12) 

The error term is uniform in yo,To, and in cq = fo/||?^o||- 

The index "I" or "11" refers to first and second order in W, but note 
that each of the corresponding contributions has an expansion in ||uo||^^. 

Proof. We first deal with the case where ||yo|| ^ 1/2- As in the proof of 
Proposition IA.21 one can integrate the equation of motion to obtain 

/■Tout ^ /"Tout 

AK {vo, yo, To) = -c / V {t)-VW {y (r) , r) dr = c / drW {y (r) , r) dr. 

J Tin J Tin 

Prom ()A.4p one easily finds, for r e [Tin,Tout], that 

||y(r)-i;o|| < c5max|r-ro| "I 

l|y(T") - (2/o + ^'o(^-'^o)) II < cc/max (r - To)^ > (A. 13) 

y{r) = y7(r) + 0(||^;o||-4) J 

where we used (jA.3p in the last hne and where 
yi (•^) = yo + vo{t - To) - c / dr / dr" VW {yo + vo (r" - ro) , r") . 

J TO J TO 

Hence 

AK {vo,yo,ro) = c T"" d^W {yi {t) ,t) dT + O {\\vo\\-'') . 

J Tin 

Expanding drW {yj (r) , r) around yo + vo{t — ro), the result follows. The 
case IIj/oII > 1/2 is handled as in the proof of Proposition IA.2i □ 
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For the purpose of proving Theorem 14. H we now turn to the computation 
of the average energy change of all trajectories with a given, sufficiently 
high, incoming momentum or energy, and that enter the ball of radius 1/2 
centered at the origin. Recalling that g{y,T) = Mg (^M~^y,iOT -\- cj)^, we 
have, for vo G K*^, b ■ vq = and k = (6, M, cf) + uitq, c), 

AE{vo,k) = AK (^vo,b-^eo,To^ . (A.14) 

We first compute the average of AE {vq, b, M, (p, c) over (p: 

/ d(l) AE{vo,b,M,(l),c) . (A.15) 

Proposition A. 4. Suppose HypothesisU\is satisfied. Then, for all vq G M*^ 

and for all 6 G M^, 6 • = 0, M G SO (d, M) , c G [-1,1], 

d4> {AE{vo,b,M,^,c) + AE{-vo,b,M,^,c)) = 

2/?}^)(eo,6,M,c 



+ 0(||^;o|r^) .(A.16) 



Here 



p\'f{eo,b,M,c) = ^ i dcp I dX I dX' {X-Xf 



T™ Jo 



xdrS (^X - ^, M, drS (^X' - 1, M, </)^ , (A.17) 

with S in, M, c/)) = VW (M-i (6 + /xeo) , • 

We remark that AE {vQ,b, M,(j),c) and AE {—vo,b, M,(j),c) are the en- 
ergy changes undergone by two distinct particles, both impinging at the 
same time on the same obstacle with the same impact parameter, but with 
opposite velocities. According to Proposition IA.31 each of the two terms 
AE {±vo,b,M,cl),c) is of order ||fo||~"'^ so that Proposition IA.4I shows that 
combining a time average with a "time reversal" vq —vq diminishes the 
energy change undergone by the particle in a scattering event drastically. 

Proof Using (|A.14p . we write, as in (|A.lip - (|A.12p 

AE = AEi + AEii + O {\\vo\\-^) . (A.18) 

It is then immediately clear from (jA.lip that 

/ d(j) AEi{vo,b,M,(j),c) = -r^ [ dX [ d(j) 
Jt™ IpoII J J 

xdrW (b+(^X- eo) , u;to + ^ 0, 
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since dr = uj ■ V,^ and W is ^-periodic. 
We now turn to AEjj and write (see ([3 



AEn = Pj^ + Pjr^+0 {Wvoin . (A.19) 

Ipoll Ihoir 



One readily finds 



p'jf {eo,b,M,(l),c) = -c^ dX dX' dX" 

Jo Jo Jo 

drS (^A - ^, M, 0^ • 5 (^X" - ^, M, 0^ . (A.20) 
and, immediately performing the i;^-average, 



^ [ f'dxf'dX' {x-xf 

^ JT™ JO JO 

drSiX- l,M, A drS fx' -l,M, 



2' ' V 2 



Note tliat in ()A.20p . the integrand is in general no longer a gradient in the 
variables, except in the special case where W (y, (p) = w (y) f ((/>). So there 

(3) 

is no reason why the (/>-average of (eo, b, M, (p, c) should vanish. But now 
remark, using (|A.20p and the definition of 5, that 



When performing in this last expression the succession of changes of variable 
defined by ^ — /i" = A — ^, = 1 — A', —fi + 1 = A", one finds 

p'ff {-eo,b,M,(l),c) = -c^ [ dX" [ dX' [ dX 

Jo J\" J\' 



drS (^X" - i, M, • 5 (^A - ^, M, . (A.21) 



Note that the domain of integration is the same as in (IA.20p . just the order of 
integration is different. So adding (|A.20p and ()A.2ip the integrand becomes 

s(x-^, M, A drS (x" - ^, M, A +drS (x - i M, A S (x" -^,M,<p\ 
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which is a total time derivative. Computing the ^-average of the sum there- 
fore yields 

d0 (eo, b, M, 0, c) + (3fj^ i-eo, b, M, 0, c)) = 0. (A.22) 
A similar computation shows that 

del) ^'ff (-eo, b, M, 0, c) = / d0 (eo, 6, M, c) . (A.23) 



Adding the various contributions, the proposition now follows from (jA.lSp . 

□ 

We are now ready to prove Theorem 14.11 
Proof of Theorem\4l\ (i) Noting that VW {y,T) = MVW {M-^y,uT + (p) 
one finds 



+ 00 



dA / d6 VVF I I A- - I e,r I = 

— oo J e-b=0 



1 



since W has compact support in its first variable. So a(^) = OSimilarly, inte- 
gration over the (p variable leads to the vanishing of a(^) (e). To prove (14. 2p . 
we first point out that, in view of the rotational invariance of the system, 
AE{M'vo,M' b,M'M, (p, c ) = AE {v o,b,M,(t),c) , for ah M' e SO(d,R). 
Consequently, AE {M'vq) = AE (vq), where ~ is defined in (j3.12p . As a re- 
sult, AE (vq) depends only on ||fo|| and not on eo. In particular AE (—vq) = 
AE (vq). It therefore follows from (IA.16P that 



A^(^;o) = ^ + O(||t;o||-^). 

This proves the first equation in ()4.2h . The second is obtained similarly and 
equation (j4.5p then follows immediately. It remains to show (j4.3p and (j4.4p . 
For that purpose, we compute B: 

^ JT'" J§<* Jb-eo=0 Jo JO 

s(b+(x-l] coa] -sfb+lx'-l) eo. 



2J V V V 2 
where S {y, cj)) = VdrW {y, 4>) . Using ([04]l below, we get 

B=^i d<j) [ dy I dy' \\y - y'\\^~'^ S {y,(t)) ■ S {y' ,(t)) , 

from which, using the definition of S and integrating by parts, we conclude 

B = id-3)^ [ d(t> [ dy f dy'\\y-yf-''drW{y,<i))drW{y\(f). 
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Remarking now that 

o2 



xdrW (^b + (^y - ^) eo) O {\\vo\\-^) 

and using again (|A.24p we finally find (|4.3p and (j4.4|) 

(ii) Note that when Vn is defined by (13. ip . k„ is independent of f„ since 
the latter only depends on for k < n. It follows therefore from ()4.5p that 

^/3n^''^ = = (^Pn^'^ — B. The same remark applies to the computation 
of the correlations. For example, when computing ((3n^(3^^f}) , for some 



positive k, one can integrate first over which yields the result because 

of (i). □ 

Lemma A. 5. In dimension d > 2, for all f : M'^ X M<^ X M+ ^ M, such that 
ho - yX^^fiyo,yo^ Wvo - y'oW) e (K''^), we have 

/g, dn (eo) /,.,^=o d6 4 dA dX'f {b + Xeo, b + X'eo, \X-X'\) 

= /iRd dyo J^a dy'o \\yo - yof "'^ / (yo, y'o, \\yo - y'oW) ■ (A.24) 

Proof. Let yo and t/q be in M'^, yo / ^o- Then there exists a unique 
(A, A', eo, 6) G M X M X S"' X M'^ with eo • eo = 1 such that 

yo = b + Xeo, y'o = b + X'eo and 6 • eo = 0. 

Since eo G S'^, there exists unique angles (^i,--- ,6*^-1) G [0, vrj'^^^ x [0,27r] 
such that 

eo = Rui, R = Rd-i {9d-i) ■ ■ ■ Ri {9i) , 

where {u\, . . . , u^) is the canonical basis of and Ri {0) is the rotation of 
angle 6 in the plan defined by ui and iXj+i. Since is orthogonal to ui, 

there exists also a unique (^p,62,-'' ^ x [0,7r]^^'^ x [0,27r] such 

that 

b = pRRu2, R = Rd-i {~ed-i] ■■■R2 {~e2, 



This gives the following equality: 

d-l d-l 

Ayody'o = \J\ dAdA'dp J] dOi \[ de„ 

i=l j=2 



where 

Rui Odxi RRu2 N M 
Orfxi Rui RRu2 N M' 



\J\ 
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with N = pRVg \^Ru2 

M = Ve(^R (j)Ru2 + Aui)) and M' = Vg (^R (^pRu2 + A' 



Ui 



Simple manipulations on the rows and columns yield 
1^1 



Rui RRu2 N 




Odxd 


Rui M' - M 



with Ji = \Rui;V eRui\ , J2 
noticing that 



ui; Ru2; V nRu2 



, the result follows upon 



d-l d-l 

dVt (eo) = ^1 n = p'^^'^.hdp JJ dOj. 

i=2 



□ 

Proof of Theorem \5.1[ Computing R (v, k) to second order in perturbation 
theory, one finds 

R{v,k) = Ri {v,k) + Rij {v,k) , (A. 25) 



where 



Ri {v, k) = -—f- / dp g \b + ne,To + 



and 



(3. f 

Rii {v, k) = j dp [K (e, K, p)-V]g{b + ^^e, tq) + O {\\v\\~^) 



with 



K (e, K,p) = / d^' / dp"g {b + p"e, tq 



—00 J —00 



As a result of ()5.ip . Rj {v, k) = 0, immediately implying a(^) = for £ = 1,2 
and hence /3W = (see = 0,1. To compute we need e ■ a^^\ 

From (|X25]) we find e • a^^^ (e, k) = T (e, k) with 

d/x y y [5(^ + /^"e'^o) • V] (e-5)(& + ;ue,ro) 

d/x / d/i" (/i-A^") [5(6 + /i"e,ro) • V] (e-5)(6 + /xe,ro) 



Noting that the integrand is unchanged under the change of variable e 
—e,p = —p,p" = —p", one finds 



[ dn (e) / dbT (e, k) 

J Jb-e=0 



' / dn (e) [ db [ dp [ dp" [p - p") 

J Jb-e=0 J -I Jjl 

X [5 (6 + p"e, To) • V] (e • g) (6 + /2e, tq) . 
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Adding the last two formulas and using the change of variables formula 
(|A.24p . we conclude 

e • a*^^) = ^'^ j '^y^y"\\y-y"\\^''^ Yliy~ y")j (^"' ^) ■ ^] {y, 0) , 

so that a partial integration yields 

= dcP [dydy"\\y-y"r''' 

X {{y-y") ■9{y,(P)) {{y-y") ■ 9{y",4>)) 
--^ [ <i<P [ (^yiy'h -y"\\^~'^9{y,(f>) ■ 9 {y",(t>) ■ 

From ()3.4p one easily sees the second term equals — ^a(^) • a^^) (e) so that 
we find, using (jS.Sp 

Wie) = [ [dydy"\\y-y"\\-'~' 

^{{y-y")-9{yA)){{y-y")-9{y"A)) 

□ 
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